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We investigate the threshold-enhanced QCD corrections to the cross sections for direct top quark 
productions induced by model-independent flavor changing neutral current couplings at hadron 
colliders. We use the soft-coUinear effective theory to describe the incoming massless partons and 
use the heavy quark effective theory to treat the top quark. Then we construct the flavor changing 
operator based on the above effective theories, and resum the large logarithms near threshold arising 
from soft gluon emission. Our results show that the resummed QCD corrections further enhance the 
next-to-leading order cross sections significantly. Moreover, the resummation effects vastly reduce 
the dependence of the cross sections on the renormalization and factorization scales, especially in 
cases where the next-to-leading order results behave worse than the leading order results. Our 
results are more sensitive to the new physics effects. If signals of direct top quark production are 
found in future experiments, it is more appropriate to use our results as the theoretical inputs for 
extracting the anomalous couplings. 
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The production and decay of top quark through FCNC 
couplings are very sensitive to new physics effects, which 
have been extensively studied in some new physics mod- 
els in the literature (see Ref. jJl and references there in) . 
Since we do not know which type of new physics will 
be responsible for a future deviation from the SM pre- 
dictions, it is necessary to study the top quark FCNC 
processes in a model independent way by an effective La- 
grangian. For t-q-g anomalous couplings, the direct top 
quark production is the most sensitive process. The anal- 
ysis based on the leading order (LO) cross sections [3 
suggests that the anomalous couplings can be detected 
down to 0.019 TcV^^ for q ^ u and 0.062 TeV^^ for 
q = c at the Tevatron Run 2. Studies with a fast detec- 
tor simulation for ATLAS indicate a similar reach at the 
LHC 0. 

As we know, the LO cross sections for processes at 
hadron colliders suffer from large uncertainties due to 
the arbitrary choice of the renormalization scale (/i^) and 
factorization scale (/.*/), and are not sufficient for the ex- 
traction of the anomalous couplings from experiments. In 
general, a next-to-leading order (NLO) QCD calculation 
is capable to reduce the scale dependence significantly. 
But it is not the case for direct top quark production 
considered here. In Ref. 0, the cross sections for the 
direct top production are calculated at NLO in QCD. 
Their results showed that the NLO QCD corrections re- 
duced the scale dependence of the total cross sections 
at the Tevatron Run 2. However, the scale dependence 
was not improved at the LHC, and even became worse 
in the region ^i^ = < mt- It was pointed out Q that 
the scale dependence comes mainly from the terms pro- 
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portional to (5(1 — z) in the partonic cross sections, and 
higher order effects are necessary to further reduce such 
scale dependence in the case of /i^ = /i/ at the LHC. 
In this paper, we report that the threshold resummation 
effects can remarkably improve the scale dependence of 
the cross sections. 

We consider the process A+B — > g{ki)+q{k2) t+X, 
where A, B are the colliding hadrons. At the parton level, 
the process is induced by the eS^ective Lagrangian 

C,^ = -9s E ^i^''''T%f!^ + ih<i^^,)qGl, + h.c., (1) 

q—u,c 

where A is the new physics scale, are the Gell-Mann 
matrices satisfying Tr(r°T'') = 6°'^/2, G^^ are the field 
strength tensors of gluon fields, n is normalized to be real 
and positive and /, h to he complex numbers satisfying 
\f? + \h? = l. 

The NLO results contain terms like (^3^) and 
■p-^jj^, which arc singular near the kinematical thresh- 
old z ^ 1. Here z = / s, = and s — {ki + fc2)^- 
Physically, these singular terms originate from the emis- 
sion of soft or coUinear gluons. The soft-coUinear effective 
theory (SCET) i5| is a natural framework to deal with 
the physics of soft and coUinear gluons and quarks. In the 
following, we will use our NLO results to derive a thresh- 
old resummation formula for direct top quark production 
using the SCET. 

We introduce the parameter ~ 1 — z 3> Aqcd/Q 
in SCETi and match the full theory operator iQJ onto 
the operator in the effective theory at leading order in 
A. This step is similar as the ones for the deep inelastic 
scattering (DIS) process in Ref. 6] and the Drell-Yan 
process in Ref. J}. However, there is a non-trivial point 
in our problem concerning the top quark field. In general, 
the top quark field can not be described in the SCET. 
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However, in the threshold region z ^ 1, the top quark is 
nearly at rest, which implies that its momentum can be 
written as q — mtv + fc, where the residue momentum 
k ^ QX^ comes from its interactions with soft gluons. 
This situation is much similar as the bottom quark in B 
decays. Therefore, we can treat the top quark field using 
the heavy quark effective theory (HQET) by a systematic 
expansion in l/mt- Thus the effective operator can be 
written as 



-39 ~ gs-^tyT'^Bnfj.Cgq{Q ,/i)W^ft^n 



'A 



(2) 



where F'^ = \ [f + ih"f^)nv(j'^'' ^ ty is the heavy quark field 
in HQET describing the top quark, is the collinear 
quark field, B!^ is related to the field strength tensor of 
the collinear gluon field in SCET and Wn denotes a Wil- 
son line which are required to ensure gauge invariance of 
the operator, n and n are two light-cone vectors satisfy- 
ing — = and n-n = 2. Cgq is the matching coeffi- 
cient which are obviously equal to unity at tree level. We 
can determine the 0{ag) matching conditions by evaluat- 
ing the on-shell matrix elements of operators in both full 
theory and SCET. In dimensional regularization (DR), 
the facts that the IR structure of the full theory and the 
effective theory is identical and the on-shell integrals are 
scaleless and vanish in SCET and HQET imply that the 
IR divergence of the full theory is just the negative of the 
UV divergence of SCET. Thus, from the 0{as) virtual 
corrections, we can obtain the NLO matching coefficient 
and the anomalous dimension of the effective operator 
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where /3o = 11/4 — n//6. 

Next, we perform the usual field redefinition to decou- 
ple the collinear and usoft fields and match the operator 
to SCETii. This matching is trivial: the matching coef- 
ficient is unity and the anomalous dimension is the same 
as Eq. and we will still denote the new operator as 
Tgq. Now we calculate the cross section in SCETn with 
Tgq{iJL) at the scale /i ~ QA^, and then match the result 
onto the product of two PDFs, which means 



CTSCET = M{z,^l) ® [fg{z,fi) ® fq{z,fl)] 



(5) 



where /iA{z,fi) is the matching coefficient. Since the vir- 
tual corrections in SCET vanish at NLO in DR, ctscet 
can be obtained by calculating the real gluon emission 



diagrams in SCET and taking into account the renor- 
malization of the operator. The result is 
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where (Tq = {Stt'^ /3)as{K/ A)'^z. Matching the cross sec- 
tion onto the product of two PDFs, all the poles are 
cancelled by the renormalization of the PDFs, 



(T 1 
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The resulting matching coeffcient is given by 
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To factorize the convolution in Eq. lO, we make Mellin 
transformation into moment space and neglect terms that 
vanish in the limit N ^ oo. The moment of A4 is 
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where N — Ne'"^ and 7_b is the Euler constant. To avoid 
large logs we choose the matching scale ji ^ Q/N. The 
running of the moments of T is governed by it anomalous 
dimension 
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Combining the above results, we can write down the 
resummed cross section in moment space 



^SCET 



(70 (m) |Csg((3^M)^ [1 + A^Ar(/^)]-^Jv(/i) 

O-o(Mr) \Cgq{Q'^,^ir)\^ 

X [1 + Mn{Q/N)\ e'^TNifif), (10) 



with 
h 



Q/N M 



^2,,(,), 



h = 



— 72(m), (11) 



where fi,. and fj,f correspond to the renomalization and 
factorization scales in the full theory, respectively, which 
are expected to satisfy /i^ ~ /■*/ ~ Q ^ Q/^ ^ Aqcd- 
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To reach the accuracy of next-to-leading logarithms 
(NLL) in the exponents, we need the coefficient of 
ln(/z^/Q^) in 71 (^) and that of In TV in 72 (m) up to two- 
loop. It can be shown that 71 and 72 have the form |^ 



71 (^) = Ai{as)ln.-^ + Ao{as), 



72(Ai) =Bi(a,)ln7V + Bo(a,), 



(12) 



and 4Ai{as) = —Bi{as)- Thus we can extract Ai from 
the information of Bi at the two-loop level. Bi is related 
to the coefficients of 1/(1 — z)^ in the DGLAP split- 
ting kernels. Using the two- loop splitting functions 
together with Eq. Q and © , we can get the coefficients 
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where Ai 




and similar for Bi^ Aq and 



J5o- Note that the coefficients Ai and Bi only depend on 
the initial states and can be related to the quark-quark 
coefficients and the gluon-gluon coefficients in a straight- 



i(i) 



forward way. For example, A\ 
{Cf + Ca)/2. 

Now we can evaluate the integrals in Eq. pi|) using 
the two-loop evolution of as in the MS scheme. Keeping 
only terms up to NLL in the exponents in Eq. H1Q(I . we 
obtain 



h+h^ lniVg(i)(/3oa,(M.)lniV/^) 
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which are the same as the results obtained within full 
QCD H [13 except the coefficients A and B. The NLL 
cross section in moment space is then given by 



^NLL 
'N 



^ aQC{as,Q'^,fir,f^f) exp (g^^) In -f g^^^) J^Nif^f), 



where C represents the contributions from Cgq , Mn, as 
well as the terms neglected in the exponents. To obtain 



the physical cross section, we perform the inverse Mellin 
transformation back to the x-space 
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We use the tricks introduced in Ref. to evaluate the 
A^-integral numerically. Finally, The resummed cross sec- 
tion at NLL accuracy is defined to be the NLL cross sec- 
tion plus the remaining terms in the NLO result which 
are not resummed, i.e.. 



„Rosum „NLL 



^NLO ^NLL 

-a —(J 



da 



NLL 



da. 



We now present the numerical results of the NLL 
threshold resummation for the direct top quark produc- 
tion at the Tevatron Run 2 and the LHC, respectively. In 
our numerical calculations, we take the top quark mass 
nit = 174.3 GeV 12], which is different from the value 
178.0 GeV used in Ref. 0. We use both the CTEQ6 
PDFs [3 and the MRST2004 PDFs 3 to evaluate the 
total cross sections. 

In Table m we list the LO, NLO and resummed cross 
sections for the renormalization and factorization scales 
fJ-r = fJ-f = nT-t using the two PDF sets. It can be seen 
that the threshold resummation effects further enhance 
the NLO cross sections by about 30%. When the re- 
summed results are compared with the LO ones, the en- 
hancement can even reach 100% in most of cases. We also 
note that the discrepancies between the different PDF 
sets are somewhat larger in the case of resummed than 
the case of LO or NLO, especially at the LHC. This may 
be due to the different fitting methods used in the PDF 
sets. We believe that when the LHC data be used in 
the global fitting of the PDFs, these discrepancies will 
be improved. 

Despite the uncertainties from the PDFs, we can es- 
timate the theoretical uncertainties by investigating the 
sensitivity of the results to the renomalization and fac- 
torization scales. We define R as the ratio of the cross 
sections (LO, NLO, resummed) to their values at central 
scale, fir = fJ-f = "Tit, always assuming fir = fif = fi for 
simplicity. Figure ^ and Figure |2 show the ratio R as 
functions of the renormalization and factorization scales 
for subprocesses gu ^ t and gc ^ t at the Tevatron Run 
2 using the two PDF sets, respectively. We can ffiid that, 
for both subprocesses, the threshold resummation effects 
further decrease the scale dependence of the NLO cross 
sections remarkably. For example, the variations of the 
cross sections in the region mt/2 < fi < 2mt, for subpro- 
cess gu t, are about 28% for NLO cross section and 
only 9% for resummed one. In the case of gc — > t, they 
are about 11% and 2%, respectively. 

In Figure O and Figure 01 the ratio R is shown for the 
cases at the LHC. First, we note again that the NLO 
corrections do not reduce the scale dependence of the 
LO cross sections. In fact, the variations of the LO and 
NLO cross sections between mt/2 and 2mt are nearly the 



4 



subprocess 


PDF 


LIIC 


f 'A A 


pb 


ToA'atroii 




•2 

fl) 


^ D.Ui LcV 1 J 


^ U.OiTcA' i J 


LO 


NLO Resum 


LO NLO Resum 


gu->-t 


CTEQ 
MRST 


12.9 
12.2 


17.0 
16.3 


23.7 
19.5 


268 425 547 
262 426 520 


gc^t 


CTEQ 
MRST 


1.71 
1.68 


2.53 
2.38 


3.71 
2.92 


13.1 28.1 38.2 
17.0 30.3 38.6 



TABLE I: The LO, NLO and resummed cross sections for direct top quark production at the LHC and Tevatron Run 2. Here 





FIG. 1: The ratio R as functions of the ronormalization and 
factorization scales at the Tevatron Run 2 for subprocesses 
gu ->■ t (left) and gc t (right) using CTEQ6 PDFs. 



FIG. 4: The ratio R as functions of the rcnormalization and 
factorization scales at the LHC for subprocesses gu ^ t (left) 
and gc -» t (right) using MRST2004 PDFs. 




FIG. 2: The ratio R as functions of the renormalization and 
factorization scales at the Tevatron Run 2 for subprocesses 
gu t (left) and gc t (right) using MRST2004 PDFs. 




It/mt iJ,lmt 

FIG. 3: The ratio R as functions of the rcnormalization and 
factorization scales at the LHC for subprocesses gu t (left) 
and gc t (right) using CTEQ6 PDFs. 



same: 18% for gu ^ t and 12% for gc t. After tak- 
ing into account the threshold resummation effects, the 
scale dependences are significantly reduced: the above 
variations become 5% and 6% for gu ^ t and gc — > t, 
respectively. Especially, the resummed results have good 
control over the scale dependences in the region n < rrit, 
where the NLO results behave worse than the LO ones. 

In summary, wc have calculated the NLL threshold re- 
summation effects in the direct top quark productions in- 
duced by model-independent FCNC couplings at hadron 
colliders in the framework of SCET and HQET. Our re- 
sults show that the resummation effects increase the total 
cross sections by about 30% in general compared to the 
NLO results. Moreover, the resummation effects signif- 
icantly reduce the dependence of the cross sections on 
the renormalization and factorization scales. Especially, 
the resummation effects have good control over the scale 
dependences in the region fj, < mt at the LHC, where the 
NLO results behave worse than the LO results. If future 
experiments at hadron colliders find signals of the direct 
top quark production, the anomalous couplings can be 
extracted by measuring the cross sections and compar- 
ing with the theoretical predictions. Thus, the precision 
of the extracted anomalous couplings directly depends on 
the accuracy of the theoretical predictions. Since the re- 
summation effects increase the cross sections and reduce 
the theoretical uncertainties, our results are more sensi- 
tive to the new physics effects and it is more appropriate 
to use our results as the theoretical inputs. 
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